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We consider the calibration of an optical quantum gyroscope by modeling two Sagnac interferom¬ 
eters, mounted approximately at right angles to each other. Reliable operation requires that we 
know the angle between the interferometers with high precision, and we show that a procedure akin 
to multi-position testing in inertial navigation systems can be generalized to the case of quantum 
interferometry. We find that while entanglement is a key resource within an individual Sagnac inter¬ 
ferometer, its presence between the interferometers is a far more complicated story. The optimum 
level of entanglement depends strongly on the sought parameter values, and small but significant im¬ 
provements may be gained from choosing states with the optimal amount of entanglement between 
the interferometers. 
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I. INTRODUCTION 


Quantum metrology and quantum parameter estimation 
offer great potential improvements in precision measure¬ 
ment. Recent experiments have demonstrated quantum 
improvements in measuring protein concentration [1], 
tracking lipid granules in yeast cells [2], and searching 
for gravitational waves [3]. In optical systems, the stan¬ 
dard way to frame problems in quantum metrology is as 
a measurement of the phase of an optical signal. The 
aim is to improve the precision of such measurements 
from the classical shot noise limit (SNL) to the quantum 
mechanical Heisenberg limit (HL) [4]. It has been rec¬ 
ognized that any practical implementation of quantum 
metrology requires methods to deal with effects due to 
environmental noise and dissipation [5]. Quantum er¬ 
ror correction has been proposed to combat the effect of 
noise [6-8], and loss-tolerant metrology protocols have 
been designed and implemented to address some of the 
negative effects of dissipation [9-12]. It has been shown 
that the measurement of d phases in an interferometer 
can obtain an improvement of a factor 0{d) in the pre¬ 
cision when multi-mode entanglement is used [13]. This 
behavior persists in the presence of photon loss [14], even 
though multi-mode entanglement is highly susceptible to 
such processes [15]. When the loss parameters are also 
estimated, there is a trade-off between the attainable pre¬ 
cision of the phase estimation and the estimation of these 
parameters [16]. However, loss and noise are not the 
only causes for imperfect metrology. The accuracy of 
a composite sensor system is only partially determined 
by the precision of the individual measurements. Other 
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sources of imperfection can include badly characterized 
responses to non-standard stimuli, or couplings between 
the parameters of interest. The performance of any larger 
scale system—i.e., one containing a number of individual 
sensors—will be limited by the presence of such nuisance 
parameters, but this aspect of quantum metrology has 
been somewhat overlooked. 

In this paper, we address the problem of nuisance pa¬ 
rameters arising from unwanted couplings between sen¬ 
sors in practical quantum parameter estimation. Such 
couplings affect the measurement precision—defined by 
the mean square error (MSE)—and must also be esti¬ 
mated, even if we are ultimately not interested in their 
numerical value [17]. For a single parameter, the quan¬ 
tum Cramer-Rao bound (QCRB) puts a lower limit on 
the MSE, determined by the inverse of the quantum 
Fisher information (QFI) [18, 19]. Multiple parameters 
lead to a QFI matrix, the inverse of which provides lower 
bounds for the MSE covariance matrix [20, 21]. Nui¬ 
sance parameters are part of this multi-parameter esti¬ 
mation problem. While the QCRB for a single parameter 
can generally be attained, this is not always true of the 
QCRB for multiple parameters [22, 23]. Where multiple 
parameters are being estimated, it matters whether the 
generators of translation of the parameters commute or 
not, with implications for the optimal strategies of the 
parameter estimation procedures [24-29]. Even though 
multi-mode entanglement can be used to improve the es¬ 
timation of multiple phase parameters beyond the classi¬ 
cal SNL [13], this is not always the case. For the example 
considered in this paper, we show that the optimum en¬ 
tanglement is a function of the nuisance parameters being 
estimated and that—for a range of parameter values with 
practical relevance—the presence of entanglement can be 
detrimental to the estimation process. 

We consider a simplified optical gyroscope configu¬ 
ration based on two Sagnac interferometers, shown in 
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FIG. 1: Two nearly orthogonal (coupled) Sagnac interferom¬ 
eters that are misaligned by an angle 6. The entire system 
rotates with angular velocity uj, and the resulting phase shifts 
ipy and ipz in the interferometers can be used to estimate oj. 
For clarity, the third interferometers measuring and the 
photodetectors are omitted. 


Fig. 1, in which the two (nominally) orthogonal inter¬ 
ferometers are misaligned by a small angle 6. We find 
that in such circumstances entanglement can hinder the 
calibration of the misalignment. In fact, we will show 
that entanglement (or classical correlations) in the quan¬ 
tum state shared by the two interferometers can limit 
the precision of the estimation process for 0, whilst en¬ 
tanglement can assist in the determination of the phases 
of interest, ipy and ipz- These results are valid for a wide 
range of physically relevant parameter values, and we in¬ 
dicate how the calibration process can be generalized to 
a set of three Sagnac interferometers measuring arbitrary 
three dimensional rotation rates. 

The misalignment of gyroscopic sensors is a well-known 
problem in the construction of inertial navigation sys¬ 
tems [30, 31]. Fibre-optical Sagnac interferometers are 
often used in modern “strapdown” inertial navigation 
systems (i.e. fixed sensors within the body of the naviga¬ 
tion system) [32]. In such systems, three fibre-optical gy¬ 
roscopes (FOGs) and three accelerometers are mounted 
in the inertial measurement unit of the navigation sys¬ 
tem. The gyroscopes provide measurements of the ro¬ 
tation rates about their axes, where the axes are nor¬ 
mally designed to form an orthogonal triad. Integrat¬ 
ing the rotation rates provides estimates of the angles 
of rotation of the system, relative to a set of reference 
axes. The angles are used to determine the system’s 
orientation, but they are also used to resolve the mea¬ 
sured accelerations into the reference axes to determine 
the system’s translational motion (i.e. velocity and po¬ 
sition). As a result, the accuracy of the gyroscopic sen¬ 
sors is often a limiting factor in the overall performance 
of an inertial navigation system. The gyroscopes will 
have mis-alignment errors due to mechanical tolerances 
in their construction and systematic errors in the mea¬ 
surement devices, both of which limit the accuracy of 
the sensors. Calibrating these errors, and correcting for 
them in software, is one way to improve the accuracy 


of the inertial navigation system, and this has become 
standard practice in many applications [30, 31]. After 
production, an inertial measurement unit will undergo 
a “multi-position” test. It is rotated through a set of 
known rotations, using a very accurate reference system, 
to obtain a static measurement value and then subjecting 
the unit to a known rotation rate after each rotation— 
normally, at least six different rotations/orientations are 
used to calibrate non-orthogonality within the triads of 
sensors, static bias measurement errors in each sensor, 
and scaling errors in the measurement of the known ro¬ 
tation rates [30, 31]. In this paper, we are primarily in¬ 
terested in the example where there is coupling between 
two non-orthogonal gyroscopes, measuring rotation rates 
about the y- and z-axes, so we will consider the simplest 
of these calibration processes, the measurement of a fixed 
(but otherwise unknown) rotation rate, followed by an¬ 
other measurement after a rotation by 7r/2 about the 
cc-axis—although we will also indicate how this may be 
extended to deal with a triad of three gyroscopes. 

II. COUPLED SAGNAC INTERFEROMETERS 

The Sagnac interferometer [33] can be described quan¬ 
tum mechanically in a very similar way to the Mach- 
Zehnder interferometer, but instead of two spatially dif¬ 
ferent paths in the latter, the Sagnac interferometer has 
a single loop with two counter-propagating modes, a and 
b. The phase shift induced by a rotation of the interfer¬ 
ometer can be written as a unitary transformation 

U (t) = exp [-iuj ■ e {ha- m) t] , (1) 

where a; is a normalized rotation rate, e is the nor¬ 
mal vector to the plane of the interferometer, and 
ha, hb are the number operators in modes a and b. 
(The Sagnac phase shift is dependent on a number 
of device-specific parameters—including operating wave¬ 
length, path length and enclosed area [32] — and it is 
proportional to the angular velocity vector applied to the 
interferometer so we will use a normalized rotation rate 
to remove to the explicit dependence on these parameters 
and to simplify the presentation of the results below). For 
simplicity, we assume that the Sagnac interferometer lies 
entirely in the xy-plane (6 = 6^), and we define ipj = uijt 
and h = ha — hb- Then we can write the transforma¬ 
tion in Eq. (1) as U{ipz) = exp[—i(pzh]. In the usual 
notation where U = eKp{—iHt/H), the Hamiltonian be¬ 
comes H = hujzh with t the interaction time (assumed 
to be known with arbitrary precision), and we will now 
set h = 1. Clearly, measurements of the phase ipz can 
be used to determine the rotation rate ujz applied to the 
gyroscope. We can construct a second Sagnac interfer¬ 
ometer in the ccz-plane (e = e^) to determine the rotation 
rate ujy, and a third can be added to determine uJx- The 
use of three such gyroscopes allows a general rotation rate 
about an arbitrary axis to be determined [34]. We con¬ 
centrate on the case with two interferometers for clarity. 
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FIG. 2: The coordinate system and the normal vectors to the 
Sagnac interferometers. Gyroscopes O and @ take one set 
of measurements, and after a rotation of 7 r /2 + S, gyroscopes 
© and O take another set of measnrements. Gyroscope O 
rotates to ©, and © rotates to O. 

but the generalisation to three interferometers will also 
be discussed below. 

In any practical construction, the two Sagnac interfer¬ 
ometers will not be perfectly perpendicular (and when 
the interferometer is constructed from optical fibres it 
may not lie perfectly in a plane). Let Uy be the number 
difference operator for the counter-propagating modes of 
the interferometer in the xz-plane, and the equivalent 
operator for the interferometer in the xy-plane. Further¬ 
more, let 0 be the angle with which the (pz interferometer 
is misaligned, shown in Fig. 1: 

§(, = cos OBz + sin 0 By . ( 2 ) 

The transformation of the optical state inside the inter¬ 
ferometers then becomes 


We will show in Sec. Ill that this leads to a linearly in¬ 
dependent set of four estimators, but first we establish 
the coordinate system local to our gyroscopes, shown in 
Fig. 2. In the original position of Fig. 1, the y-rotation 
0 is matched to the local y-axis by definition, while the 
0 -rotation O is misaligned according to Eq. (2). After a 
rotation of 7r/2 -|- (5 about the sensor a;-axis, gyroscope @ 
is now aligned along the normal e'y with 

e'y = cos{0 + d)ey — sm{0 -I- <5) e^ , (5) 

and gyroscope O is aligned along the normal e" with 

e” = — cos Sez — sin S By . (6) 

The measured phases are decomposed in the same way 
as the normal vectors {(pj -f-)- e^). In the ideal case where 
0 = (5 = 0 we have Lpy = ip'y^ and ipz = pi'z = —^'z- 

Let the joint state in the two nearly perpendicular 
gyroscopes be denoted by pij, where i and j indicate 
the gyroscopes 0-0 in Fig. 2. The two gyroscopic 
measurements are equivalent to a single measurement 
with four gyroscopes simultaneously, with a joint state 
p = pi 2 ® P 34 . In general, the optimal state pi 2 for O 
and @ will not be the same as the optimal state P 34 for @ 
and O, since the evolution 17 i 2 of the gyroscopes O and 
© is not equal to the evolution U 34 for the gyroscopes @ 
and O due to the different relative rotation u>. The trans¬ 
formation of the optical state due to the rotation rate u} 
can be written on the joint system as U = 1/12^134, with 

17i 2 = exp [-i {(p'zfiz + p-yfiy)] (7) 

U 34 = exp [-i {ip'yfi'^ + P’'zn'y)] , (8) 

which depend in a nontrivial way on the parameters 93 . 
The four operators hy, fiy, fiz, and hi, commute, and can 
be measured simultaneously. The joint evolution then 
becomes 


U{(p) = exp [—i {(pyiiy + cos 0 (fzfiz + sm0 ipyhz)] , (3) 
leading to a Hamiltonian for the system 

H = UJyfly + COS 0 UJzflz + sin 0 UJyflz ■ (4) 


U{ip) = exp|-i ^(l)yny - /3 (^ + ^) } 

X exp j-I ^ + I) + 13(0 + (5)hi,j |(9) 


where 


There is now a coupling between the two interferome¬ 
ters given by the term sin0ujyfiz- As a consequence, we 
have three unknown parameters, ipy, ipz, and 0 , but we 
measure only two observables, fiy and fiz- The problem 
is therefore underdetermined, and we cannot extract the 
true values of py^ pz without an unknown bias. 

To remedy this, we may rotate the system of gyro¬ 
scopes by 7r/2, whilst keeping the applied rotation rate 
fixed with respect to an external reference system. This 
allows us to measure four different observables (two for 
each orientation), given our three unknown parameters. 
Assuming that the system is rigid, 0 remains unchanged. 
However, we do need to introduce a new nuisance pa¬ 
rameter 5 that encodes imperfections in the 7 r /2 rotation. 
This leaves us with four parameters and four observables. 


(}{a) = (j)y cos a — (jiz Sana . ( 10 ) 

and p> = {py,pz,0,S). The evolution U in Eq. (9) is 
expressed entirely in terms of the measurable observables 
fly, fly, fiz, and fi'^, and the four unknown parameters ip. 
In the next section we use this evolution to calculate the 
quantum Eisher information and the Cramer-Rao bound 
for these parameters. 

III. COVARIANCE AND FISHER 
INFORMATION 

To determine the ultimate precision with which we can 
estimate the Sagnac phases and the couplings between 
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them, we consider the quantum Cramer-Rao bound 

Cov(9?) > ^ , (11) 

where Cov(c^) is the covariance matrix of the four vari¬ 
ables y, N is the number of independent measurements, 
and Iq(<^) is the quantum Fisher information (QFI) ma¬ 
trix of the three variables [18, 35] with elements: 

[Iq(v)]*j = 2 didj log , (12) 

where di is the derivative with respect to (pi, and the 

derivative with respect to (pj. Let G = {Gy, G^, Gg, G^)”*" 
be the tuple of generators of translation in our four pa¬ 
rameters. Generally, a generator of translation Gq, of a 
parameter a can be defined as [36] 

G„ = (13) 

This allows us to relate the derivative of the quantum 
state j^) with respect to pi to the generator Gi via a 
Taylor expansion of U. Evaluating the matrix elements 
of the QFI matrix for pure states j'0) then yields 

IWh = 4 ({G„G,})^ - (G,)^(G,)^) 

= 4[Cs(G)].,, (14) 

where {0)cp = {ip{ip)\0\4>{‘p)) for some operator O, and 
[Cs(G)]ij is the symmetrized covariance matrix element 
between operators Gi and Gj, originating from the fact 
that the quantum Fisher information matrix in equation 
(12) is derived from the symmetric logarithmic derivative 
[18]. Since all our generators commute with each other, 
we can ignore this technical requirement and drop the 
subscript S. 

Applying Eq. (13) to U{(p) in Eq. (9) for the parame¬ 
ters py, pz, 0 and S, we obtain the generators 

Gy = fly + sin 0 hz — sin 6 fi'y + cos{0 + S) , 

Gz = cos 9 fiz — cos 5 fi'y — sin(0 -\- 5) h'^, 

Gg = P{9) fiz + 13 (^9 + 6 + —^ > 

Gs = -m n; + /? (0 + d + I) h; . (15) 

The relation between the generators Gj and the observ¬ 
ables fi]^ is linear and can be expressed in matrix form as 
G = Mix with 

/I sin0 cos(0-I-(5) — sin(5\ 

0 cos9 sin(0-I-(5) — cosi5| , . 

“ 0 /3(0) /3(0 + (5 + 7r/2) 0 ’ 

Vo 0 /3(0 + (5 + 7r/2) -(3{5)) 

and n = (h^, h^, h^, h),)”*". The determinant of this ma¬ 
trix is 

det M = [/3(0) + /3((5)] /3 (0 + 5 + 0 cos(0) 

-/3(0)/3(d)sin(0 + d), (17) 


which is nonzero for most values of p, and in particular 
for our case of interest of small values of 9 and <5. The 
four observables can therefore be used to determine the 
four parameters unambiguously. 

We calculate the QFI in terms of the matrix M : 

[lQMh.=4[(G.G,)-(G.)(G,)] 

= 4^M,feM,,[C(n)]fe,, (18) 

kl 

or 

Iq(¥ 5) =4MC(n)MT, (19) 

where [C(n)]fe/ = {fikfii) — {hk){fii) are covariances that 
depend only on the state inside the interferometers (all 
the alignment information is encoded in M). This ex¬ 
pression contains a large number of variables in C(n), 
corresponding to the extensive freedom to choose input 
states of the interferometers. However, we can drastically 
reduce the number of variables using simple symmetry ar¬ 
guments: Since the state p is a tensor product pi 2 ® P 34 , 
the 2x2 off-diagonal sub matrices of C(fi) are zero due to 
the fact that for this case {hkfii) = {hk){fii)- Moreover, 
to keep the resources in the gyroscopes identical between 
rotations, we take the photon number differences (Ahy)^ 
and (Afiz)'^ the same in the gyroscope setting 0 -|-@ and 
©-I-0. This leads to the covariance matrix 


/ (Ahy)2 

p^(12) 

C/yz 

0 

0 \ 

^(12) 

L^yz 

(Ah,)2 

0 

0 

0 

0 

{Ahyf 

^(34) 

^yz 

V 0 

0 

p.(34) 

<^yz 

{Afizfj 


where GyP = {hyfiz) — {f^y){f>'z) for gyroscope pair i and 
j. A priori there is no reason to choose different probe 
states for the two orientations 0 -|-@ and @-|-0 during 
normal operation. However, during the calibration stage 
of the gyroscopes the rotation is a precisely known 
rotation in an external reference frame, and its value will 
generally determine different optimal states pi 2 and P 34 . 
Here we choose pi 2 = P 34 to keep the analysis tractable. 
The optimal strategy for different pi 2 and P 34 , as well as 
the optimal rotation direction a; will be the subject of 
future work. 

Without prior knowledge of 9 and S there is no reason 
to require different values for (Afiy)^ and (Ah^)^. We 
can therefore take 

{Ahy)'^ = {Arizf = (An)^ (21) 

and 

Gfe) = XiAnf , ( 22 ) 

with — 1 < A < 1 the correlation coefficient. Note that 
the correlation in this context means entanglement, since 
the optimal states are pure states [19], and classical cor¬ 
relations require mixed states. 
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FIG. 3: Color online. The MSEs for ipy, ifiz, 9 and <5 as a 
function of (a) the variance in the photon number (Ah)^ in 
each Sagnac interferometer—assuming all four interferome¬ 
ters having the same (An)^ with A = 0, and (b) the corre¬ 
lation coefficient A between photon number difference in two 
Sagnac interferometers. Here, 9 = 0.02 rad, 5 — 0.013 rad, 
ifiy — 0.66 rad, ifiz = 0.17 rad, and (Ah)^ = 10. The curve 
legend for the variances of ipy, ipz, 9 and 5 applies to both 
figures. 


The diagonal elements of Cov((p) are the variances of 
the parameters of interest, namely Var</jy, Vari^z, and 
the nuisance parameters Var 9 and Var 6 . We can choose 
to optimize any one of these variances, two or three of 
them, or all four. In the latter case, we need to choose a 
quantum state that minimizes 


Tr[Cov(93)] > 


1 

N 


Tr 


ioiv) 


(23) 


The right-hand side of Eq. (23) provides a bound on the 
optimal joint estimation of ip that may be achieved in the 
asymptotic limit of large N. Another interesting optimi¬ 
sation is to minimize the combination Maxpy 
The values of the nuisance parameters 9 and 5 are only 
interesting in as far as they can be used to improve the 
accuracy of the overall sensor. They do not convey in¬ 
formation about the rotation of the gyroscopes per se, 
but—once calibrated—they can be used to correct hxed 
errors due to couplings between the measurements of the 
interferometers. The MSE for a typical case of 0 = 0.02 
rad, 5 = 0.013 rad, Py = 0.66 rad , and Pz = 0.17 rad is 
shown in Fig. 3. 


confirms our intuition that we should maximize the vari¬ 
ance (Ah)^ in each Sagnac interferometer (for example 
using NOON states [37]). 

We can also investigate the dependence of [lQ^(<^)]fcfc 
on the correlations in the states between the interferome¬ 
ters, A. Consider first a perfect gyroscope with 9 = 6 — 0. 
Solving the simpler problem of only two generators fiy 
and fiz and ignoring the nuisance parameters entirely we 
obtain 


Var py > 


Var pz > 


1 (Ah^)^ 

AN (Ah^)2(Ah,)2-C2, ’ 

1 i^ny)^ 

AN (‘Ah^)2(Ah,)2-C2, ■ 


(25) 


On the other hand, taking the limit of 0, (5 —?> 0 in the 
four-parameter problem and setting n = Pyjpz yields a 
different result for Var pz '■ 


^1 {NfiyY + K^{/S.nzY + 2nCyz 
{NhyY{AhzY - ciz 


(26) 


and this reduces to Eq. (25) only in the limit where k —> 0 
(up to an overall factor of two, since Eq. (26) contains 
two independent sets of measurements—initial and ro¬ 
tated by 7r/2—in this limit). The reason for the discrep¬ 
ancy between Eq. (25) and Eq. (26) is that the former 
does not take into account the nuisance parameters at 
all, while the latter takes into account the uncertainty in 
the nuisance parameters 9 and 5 when their mean value 
is equal to zero. 

When we consider the mean square errors in Eq. (25), 
we see that the right-hand side is minimized when 
Cyz = 0. Since these are effectively two separate 
single-parameter measurements, the quantum Cramer- 
Rao bound is attainable asymptotically, and we con¬ 
clude that, in the absence of alignment errors, any cor¬ 
relation between the input states in the two Sagnac 
interferometers—including entanglement—is detrimental 
to the precision of the optical quantum gyroscope. 


B. Correlation coefficients for misaligned 
gyroscopes 


A. Variances for well-aligned gyroscopes 


The QFI in Eq. (19) based on the generators in Eq. (15) 
and the covariance matrix in Eq. (20) is easily calculated. 
The diagonal elements of its inverse can be written as 



1 Fk{pA) 

4[(Afi)2[detM(<p)]2 (1 — A^) 


(24) 


where the functions are lengthy expressions in p and 
A (given in full in Appendix A). Note that the diago¬ 
nal elements of the inverse QFI are all proportional to 
(Ah)“^, as well as the squared determinant of M. This 


When alignment errors are introduced, the situation is 
much more complicated. Fig. 3a shows the dependence 
of the parameter errors as a function of (Ah)^. Fig. 3b 
indicates that, for the values of Py and Pz shown, the 
combination Var py F Var pz is optimal for some non-zero 
value of the correlation coefficient A, whilst the mini¬ 
mum mean squared errors in the nuisance parameters, 
9 and <5, are close to zero correlation for this example. 
This demonstrates that entanglement can be beneficial in 
some circumstances, at least for a subset of the unknown 
quantities, and the optimum correlation is not necessar¬ 
ily the same for all of the parameters. We therefore need 
to study this behavior in more detail. 
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fy 

tpz Var ipy + Var ipz 
(A = 0) 

Var ipy + Var (pz 

(A = Aopt) 

Aopt 

0.01 

0.01 

0.05 

0.0467 

-0.2679 

0.01 

-0.30 

0.0375 

0.0375 

0.0111 

0.20 

0.01 

5.0375 

5.0251 

-0.0498 

0.20 

0.20 

0.05 

0.0467 

-0.2679 

0.20 

-0.30 

0.0431 

0.0414 

0.2014 

-0.30 

0.01 

11.288 

11.275 

0.0333 

-0.30 

0.20 

0.0656 

0.0597 

0.3139 

-0.30 

-0.30 

0.05 

0.0467 

-0.2679 


TABLE I: The sum of the variances Martpy + Var</52 for se¬ 
lected values of ipy and with (Ah)^ = 10, both when 
no entanglement is present between the gyroscopes (A = 0) 
and when the optimal probe state is used (A = Aopt)- Here, 
6 — 0.02 rad, and 5 = 0.013 rad. 


When operating as a gyroscope, the input states should 
be selected to minimize the variance of the measured 
Sagnac phases, ipy and tpz, thereby improving the accu¬ 
racy of the rotation rates measured by the sensor. This is 
the asymmetric condition shown in Fig. 3b, minimizing 
the combination Var ipy + Var ipz to find the optimum cor¬ 
relation coefficient Aopt- Some example values are shown 
in Table I. The variances depend strongly on the actual 
values of Lpy and ipz but the gain in precision by choosing 
Aopt instead of A = 0 appears to be modest for these ex¬ 
amples. The optimal states for a given parameter can be 
found by constructing equal superpositions of the eigen¬ 
states with minimum and maximum eigenvalues of the 
corresponding generator of translations Gj in Eq. (15) 
[4]. It is generally difficult to create these optimal states, 
and in addition we require a strategy to choose between 
different probe states optimized for different parameters. 
Given the modest improvement in precision it is ques¬ 
tionable whether the extra effort would be merited. The 
typical reduction in the variances shown in Table I is of 
the order of 5-10%. However, it is important to remem¬ 
ber that the accuracy of gyroscopes in a strapdown iner¬ 
tial navigation affects the accuracy of all of the derived 
quantities, because the orientation information provided 
by the gyroscopes is used to resolve the measured accel¬ 
erations into a set of reference axes to then determine the 
velocity and position information. This sensitivity to er¬ 
rors in orientation can mean that even marginal gains in 
the accuracy of the gyroscopes may be important for the 
overall performance of the system. The trade-off between 
the difficulty in preparing the entangled states and the 
benefits of using these states will therefore be application 
specific. 

The case that we are considering here is not the general 
operation of a gyroscope, it is the situation where the 
gyroscope is being calibrated to estimate the misalign¬ 
ment of the individual sensors, using the multi-position 
test described above. In cases where the phases are com¬ 


pletely unknown initially, a reasonable approach is to ex¬ 
amine the behavior of this optimum value of A as a func¬ 
tion of 9 and 5 when averaged over all possible values 
for ify and ipz- These results are shown in Fig. 4a and 
Fig. 4b, minimizing the errors in the measured phases 
Var Lpy -|- Var Lpz as before. The optimum value of A for 
a wide range of the nuisance parameters is very close to 
zero (—I < 9,6 < |), when averaged. In particular, the 
optimum value is zero for small values of 9 and 6 , and 
the response of each is only very weakly dependent on the 
value of the other nuisance parameter. (A non-zero value 
for A would indicate that entanglement is beneficial in 
improving the overall precision of the optical gyroscope). 
This implies that, in the absence of information regard¬ 
ing the measured phases, the calibration of small nui¬ 
sance parameters alone will be hindered by the presence 
of correlated input states. It is only when the coupling 
between the two gyroscopes becomes significant that cor¬ 
related input states and entanglement could be beneficial. 
This is in contrast to the case where the measured phases 
are known initially. Fig. 4c and Fig. 4d, where the opti¬ 
mal value of the correlation coefficient is given as a func¬ 
tion of the nuisance parameters for (py = 0.66 rad and 
(pz = 0.17 rad. In the cases shown, the optimum value 
for A is strongly dependent on 9 and 6 . There is clearly 
a complicated relationship between the optimal degree 
of correlation between the two Sagnac interferometers, 
which varies significantly not only with 9 and 6 , but also 
with ipy and (pz- 

Starting with unknown values for the measured phases, 
the results shown in Fig. 4a and Fig. 4b indicate that 
correlated inputs will not necessarily help to improve the 
accuracy of the estimates of the nuisance parameters, for 
small values of 9 and 5 at least. However, as the es¬ 
timates of the Sagnac phases (py and ipz improve, the 
results shown in Fig. 4c and Fig. 4d indicate that the 
optimum value for A becomes strongly dependent on the 
values of the nuisance parameters. This implies that the 
optimal calibration process could involve an adaptive ap¬ 
proach, with the selected value of A being dependent on 
the estimated values and expected errors in ipy and ipz 
and the estimated nuisance parameters themselves. The 
problem becomes even more complicated when one takes 
into consideration that the optimum value of A calcu¬ 
lated in Fig. 4 is the value that minimizes the combina¬ 
tion Maripy Marifz, i.e. the variance in the phase mea¬ 
surements, not the variance in the nuisance parameters. 
Fig. 5 shows the behavior of the optimum correlation co¬ 
efficient when minimizing the variance in the estimated 
value for each of the nuisance parameters, again averaged 
over the measured phases. 

Whilst the different optimisation criteria shown in 
Fig. 5 do show qualitative and quantitative differences 
when compared with the results shown in Fig. 4a and 
Fig. 4b, the results are not inconsistent —meaning that 
the optimal values for A are zero or close to zero in each 
case, i.e. in the physically relevant regime where the 
misalignment parameters are small. It is only when 6 
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FIG. 4: Color online. The correlation coefficient A that pro¬ 
duces the lowest joint variance Var ipy + \/ar for (An)^ = 10 
as a function of a) 9 and b) 5, both averaged over all values 
for the phases, for different values of 5 and 6 respectively. 
The averaged values are calculated by averaging over values 
ol 5,9 = ±0.01, ±0.1, ±0.2 and ±0.3 rad. c) and d) give the 
same calculations for specific values of the measured phases, 
ipy = 0.66 rad and ipz = 0.17 rad, and the nuisance parame¬ 
ters. 

is relatively large (|(5| ~ 0.3 rad) that the optimal value 
for the correlation is non-zero at 9 ~ 0.0. This means 
that an adaptive calibration method should be feasible, 
although the precise form that it should take will depend 
on the criteria being optimized and the quality of the a 
priori knowledge regarding the rotation rates applied to 
the system, and the corresponding Sagnac phase mea¬ 
surements. Ideally, the calibration would be performed 
using a rotation rate that had been selected to produce a 
minimum variance for the parameter estimates, and these 
vary signihcantly, as can be seen in Table I. Fig. 6 shows 
the mean square errors in the four estimated parameters 
as functions of the applied Sagnac phases. 

C. 3D gyroscopes 

In our discussion so far we have ignored the third dimen¬ 
sion. Here we give a simple argument that shows that 
the calibration process for quantum interferometers can 
be generalized to three dimensions in a straightforward 
way. Consider the three axes determined by three Sagnac 
gyroscopes. We pick one axis as a reference axis—such 
as the y-axis in the discussion above. The remaining two 
axes deviate from a perfect triad in three angles, namely 




FIG. 5: Golor online. The correlation coefficient A that mini¬ 
mizes the variance in the estimates of the nuisance parameters 
for (Ah)^ = 10, averaged over all values for the phases: a) 
A as a function of 9 minimizing Var 6, b) A as a function of 
5 minimizing Var 5, a) A as a function of 9 minimizing Var 5, 
and a) A as a function of <5 minimizing Var 9 —each for several 
values of the other nuisance parameter and when averaged 
over values of 5, 0 = ±0.01, ±0.1, ±0.2 and ±0.3 rad. 


the angles for the x and z-axes with the y-axis, and the 
angle between the a:-axis and the z-axis. These are three 
nuisance parameters 9 (now a three-dimensional vector) 
arising from the physical tolerances inherent in any man¬ 
ufacturing process. Calibration requires that we deter¬ 
mine these angles with high precision. During operation, 
any Sagnac phase measurements can be corrected by re¬ 
moving the coupling introduced by the nuisance param¬ 
eters in post-processing. 

The calibration consists of applying a well-defined ro¬ 
tation rate to the system, fixed with respect to an exter¬ 
nal frame, which leads to expected values of the phases 
ipx, ipy, and ifz- We must estimate these phases with the 
gyroscope. We then rotate the gyroscope sensor through 
a known angle and repeat the measurements with the 
same rotation rate applied. For each independent set 
of phase measurements we take three readings, namely 
the photon number difference in x, y, and z. The ques¬ 
tion is now if we can take enough different readings by 
applying these known rotations to the three-dimensional 
gyroscope several times, taking into account that each 
finite rotation is itself associated with a nuisance param¬ 
eter, such as the parameter <5 above. If we take K sets 
of readings, after each rotation about an arbitrary an¬ 
gle, 5j (1 < j < K), we introduce—in principle—three 





































































FIG. 6: Color online. The mean sqnare errors for ip as fnnc- 
tions of a) Py with pz = 0.17 rad, and b) pz with py = 
0.66 rad, for (Afi)^ = 10, 9 = 0.02 rad, and 5 — 0.013 rad. 
The blue solid line indicates the combined log MSE in py and 
Pz, while the orange dotted line indicates the combined log 
MSE in 9 and <5. Calibration using specially selected Py and 
Pz achieve the best results when adaptive strategies are used, 
i.e. when the values of Py and pz are chosen based on the 
inferred values of 9 and 5. 


new nuisance parameters for each 8 j. It is easy to see 
that the total number of parameters stays ahead of the 
total number of readings by three, and we will not be 
able to establish all parameter values unambiguously. To 
overcome this problem, we rotate the gyroscope around 
the same externally defined axis that determine the cal¬ 
ibration signal {}Pxi ‘T’y, (^z). This will incur only a single 
nuisance parameter for every rotation (after the first set 
of measurements), namely the rotation angle, and it still 
allows us to rotate the gyroscope in any direction we wish. 
Suppose we take a set of measurements in K gyroscope 
orientations. We will have the three components of the 
applied rotation rate, three angles corresponding to the 
mis-alignment of the Sagnac axes, plus K —1 rotation an¬ 
gles, giving a total ot K + 5 parameters that need to be 
estimated from 3K measurements of the Sagnac phases. 
A successful calibration therefore requires 

3K>K+ 5, (27) 

which means that we need at least three sets of measure¬ 
ments. If we then wish to relate the gyroscope orientation 
to an external reference frame we need another three pa¬ 
rameters, and this can be achieved by taking another set 
of measurements at a fourth orientation. This analysis 
implies that the generalisation of our analysis to three 
dimensions does not pose a major obstacle—which is to 
be expected since similar processes are already used in 
inertial navigation systems and can include more compli¬ 
cated sequences of rotations and applied rotation rates to 
characterize other errors [30, 31] —and there is no reason 
to believe that the main results given above for entangled 
input states would simplify significantly in moving from 
two to three dimensions. If anything, the optimisation of 
the input states in three dimensions would be expected 
to be an even more complicated problem than the case 
presented here. 


IV. CONCLUSIONS 

We have studied the role of entanglement in the cali¬ 
bration and operation of an optical quantum gyroscope, 
relying on measurements of the Sagnac phase shift. The 
accuracy of the gyroscope is an important factor in in¬ 
ertial navigation systems, and any imperfections in the 
system must be well characterized for reliable operation 
as a rotation sensor. The interferometers are sensitive 
to rotation rates, which must be integrated with respect 
to time to generate orientation information, which can 
then be used to resolve acceleration measurements into a 
reference frame, thereby allowing changes in velocity and 
position to be determined. 

For a perfectly aligned gyroscope, the optimal state 
in each individual interferometer will be highly entan¬ 
gled, NOON states for example, with no entanglement 
or other correlations between the interferometers. The 
Sagnac phase measurements are essentially independent, 
and the quantum Cramer-Rao bound can be achieved 
asymptotically. When the interferometers are not per¬ 
fectly orthogonal to each other, the appearance of nui¬ 
sance parameters (the mis-alignment of the sensors) will 
give rise to couplings between the measured values for 
the rotation rates. Knowing the values of these nuisance 
parameters allows the effect of these couplings to be re¬ 
duced or removed, thereby improving the accuracy of the 
system. To this end, the system must be calibrated by 
taking readings from the gyroscope in a variety of known 
orientations. In two dimensions, we have demonstrated 
that the optimal states for such systems are dependent 
on the values of the nuisance parameters and the condi¬ 
tion that is being optimized. Specifically, the states that 
would be optimal for the operation of the gyroscope as 
a sensor—i.e. the states that minimize the variance of 
the phase measurements—can be different to states that 
optimize the calibration of the nuisance parameters. In 
fact, the calibration of the optical quantum gyroscope is 
complicated by the fact that the variances of our estima¬ 
tors depend strongly on the values of the various param¬ 
eters in the problem. In addition, we have found that 
the precision of the measurements of both the nuisance 
parameters and the applied signal is highly dependent 
on the amount of entanglement, as indicated by the cor¬ 
relation coefficient A. These dependencies must be fully 
understood, as the use of suboptimal probe states can 
make the calibration worse. Small improvements in the 
estimation of the nuisance parameters may seem to come 
at the expense of hard to engineer quantum probe states, 
but the nature of error accumulation in inertial naviga¬ 
tion means that even modest gains in accuracy can lead 
to an important operational improvement of slower error 
divergence. 

An important open question is how to establish the 
optimal quantum strategy for the calibration of the op¬ 
tical quantum gyroscope: what are the optimal quan¬ 
tum states to send into the three interferometers, and 
how should we divide our resources between the different 
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probe states? Even when the measured observables com¬ 
mute, it is not immediately clear how the probe states 
must be chosen in order to maximise the information gain 
about the different parameters in the gyroscope. These 
questions will be addressed in future work. 
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Appendix A: Explicit forms for lA 


Here we give the explicit forms for Fk in Eq. (24) in the case where (Ahj,)^ = (Ahz)^ = (An)^ and Cyz = A(An)^. 
For notational compactness, let 7 = /I(0 -I- J -I- 7 r/ 2 ), Pe = /3(0), and fis = The functions then become 


Fy 

F, 


Fe 


Fs 


[j{l3s + /3e) cos 9 - fisPe sin((5 -I- 6)f = (det M)^ , 

7 ^ [fij (2Asin0 -I- sin^ 0 - 1 -1) -|- 2 /? 5 ^e sin J(sin0 F X) + (sin^ J -I- l)] -|- PjjSl [cos^(J -I- 0) -|- l] 
- 2^I3sP9 {Pscos{6 + 9){sm9 + X) + (3g [sin(5cos(<5 -|- 0) -I- A]} , 

7 ^ cos^ 0 [(sin J — sin 6 *) (sin J — sin 0 — 2A) -I- 2] 



A sin(2(5 -|- 39) + sin 6 sin((5 -I- 29) + 


1 

2 


cos(2(5 -|- 49) -I- A sin 0 


(Al) 

(A2) 

(A3) 


— 7 /3^ COS0 {A [cos(2J + 9) — 3cos(J -I- 29) + cosS + cos0] -I- 3sin((5 -I- 0) -I- sin(2i5 -I- 29) — sin(i5 -I- 39) — sin 20} , 

7 ^ cos^ 9 [(sin 6 — sin 9) (sin 9 — sin S + 2A) — 2 ] (A4) 

— P 9 { (sin^ J -I- 1) sin^(J -I- 0) -|- cos^ 0 [cos^(i5 -I- 0) -I- l] — 2cos0sin(i5 -I- 0) [sin(5cos((5 -I- 0) -I- A]} 

-I- 7 cos 9^0 [A(cos(2(5 -I- 0) -I- cos((5 -I- 20) -I- cos J — 3 cos 0) -|- 3 sin((5 -I- 0) -I- cos(20 -I- 25) sin(0 — 5) — sin 25]. 
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